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AVIATION  AND  SPACE  TECHNOLOGY 


UDC  593.3 

STUDY  OF  Dm^MICS  OF  CIRCULAR  RING  STRUCTURES  ON  PROGRESS-40  CARGO 
SPACECRAFT 

917F0158A  Kiev  DOKLADY  AKADEMII  NAUK  UKR4INSKOY  SSR:  SERIYA  A  — 
FIZIKO-MATEMATICHESKIYE  I  TEKHNICHESKIYE  NAUKI  in  Russian  No  11,  1990 
(manuscript  received  14  Axig  90)  pp  36-43 

[Article  by  academician  B.  Ye.  Paton,  Yu.  P.  Semenov,  P.  M.  Belousov  et 
al.,  Kiev  Construction  Engineering  Institute] 

[Text]  A  space  experiment  to  unroll  two  large  structures  20  m  in 
diameter  each,  which  use  drives  frc»n  material  having  a  memory  effect  for 
unrolling,  was  conducted  on  the  Progress-40  spacecraft  and  the  Mir 
orbital  station  during  the  period  from  3  through  5  March  1989. 

The  process  of  unrolling  and  shaping  these  structures  under  space  flight 
conditions  was  studied  experimentally,  and  the  stiffness  and  dynamic 
characteristics  of  the  unrolled  loop  structures  during  dynamic 
operations  with  the  cargo  spacecraft  were  studied. 

The  free  vibrations  of  large  circular  structures  Eind  their  transient 
motions,  caused  by  programmed  progressive  and  rotational  maneuvers  of 
the  cargo  spacecraft,  were  calculated  in  jjreparation  of  the  experiment. 

Let  us  consider  the  problems  of  mathematical  modeling  of  the  dynamics  of 
the  carrier-ring  system,  the  center  of  mass  of  which  moves  in  a  circular 
orbit  about  the  earth.  When  performing  maneuvers  of  the  carrier 
according  to  the  given  program,  the  rings  connected  to  it  perform 
complex  motion,  consisting  of  translational  motion  of  the  body  itself 
and  the  relative  motion  of  the  rings  with  respect  to  it.  Let  us 
disregard  the  orbital  motion  caused  by  gravity  and  inertia  and  let  us 
assume  that  the  system  is  free  to  study  the  dynamics  of  relative  motion 
of  the  rings.  Let  us  introduce  the  inertial  coordinate  system  OaXYZ  and 
the  system  Ocxyz,  rigidly  bound  to  the  carrier  so  that  its  axes  coincide 
with  the  principal  central  axes  of  the  system,  and  axis  OcZ  passes 
through  the  centers  of  rod  rings,  while  axis  OcX  is  perpendicular  to 
their  plane. 

The  equations  of  dynamic  equilibrium  of  the  deformed  element  of  the  rcxl 
ring  are  represented  in  the  form  [4] 
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(1) 


dF  ,  c  ,  t  n  dM  r  n 

— — 0)  X  F  4-  f  =  0,  — - r®xM-rTXF  =  0, 

ds  ds 

vdiere  F  is  the  internal  force  vector,  M  is  the  internal  moment  vector,  s 
is  a  coordinate  measured  along  the  axis  of  the  line  of  the  rod  by  the 
distance  from  some  initial  point  to  the  current  point,  f  is  the  vector 
of  the  external  load  intensity,  viiich  also  includes  the  force  of 
inertia,  T  is  the  unit  vector  of  the  tangent  to  the  axial  line  of  the 

rod,  is  a  Darby  vector,  and  d/ds  denotes  the  local  derivative. 

Let  us  introduce  the  coordinate  system  (u,  v,  w) ,  bound  to  the 
considered  cross-section  of  the  rod,  axis  w  of  which  is  directed  along 
the  tangent  to  the  flexible  line,  and  axes  u  and  v  are  directed  along 
the  principal  central  axes  of  inertia  of  the  cross-sectional  area. 

Projections  Mu,  My,  and  of  moment  M  onto  the  corresponding  axes  are 
expressed  by  the  stiffnesses  upon  bending  of  A  and  B  and  by  twisting  of 
C,  and  also  of  curvature  p,  q  of  element  ds  and  twisting  r  of  the  rod 
axis  by  the  relations 

■^{p —  Po)’  ^v  =  B{q  —  qa),  =  C(/'  — /■„),  (2) 


in  vdiich  the  values  of  p,  q,  and  r  for  a  deformed  state  of  the  rod  are 
den toed  by  zero. 

Let  us  supplement  system  of  equations  (1)  and  (2)  kinematically  by 
Frenier  equations  [1] 


dx  __  jT_  _  JL  _i_  "  ( 3 ) 

ds  ~  p  ’  ds  p  '  y  ’  y  -  —  jT  ■ 

Here  n  and  b  eire  unit  vectors  of  the  normal  and  binormal  and  p  eind  T  are 
the  radii  of  curvature  and  twisting. 

Let  us  also  write  the  equations 


dx  dy  _ 

ds  ~  ds 


X,. 


(4) 


The  system  of  solving  equations  (l)-(4),  formulated  in  the  domain  of 
variation  0  <  s  <  S  of  independent  variable  s  has  a  common  18th  order. 
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The  first  six  integrals,  which  follow  from  the  conditions  of 
orthonormalization  of  the  Frenier  basis, 

|xl=  1,  lnl=  1,  T-n  =  0,  TXn  =  b 


close  the  system  of  solving  equations  and  permit  us  to  avoid 
redetermination  of  it. 

To  construct  the  vibration  equations  of  the  rod,  let  us  use  the 
D’Alambert  principle  and  let  us  represent  the  force  acting  on  it  in  the 
form 


f  =  —  7a, 


(5) 


where  7  is  the  linear  density  of  the  rod  and  a  is  the  acceleration 
vector  of  the  element  of  the  rod  in  system  OaXYZ. 

Let  us  use  the  Coriolis  formula  [2]  a  =  a^,  -f-  a^.  -j-  a^.,  vAiere  Sr  =  x  I  +  j/ j  -j-  2k 

is  relative  acceleration,  a,,  =  +  exp  +  ^2x(Qxp)  is  translational 

acceleration,  a^  =  2Qx(.v'i -j- +  zk)  is  Coriolis  acceleration,  ao  is  the 

acceleration  of  the  center  of  mass  of  the  carrier ,  0  and  c  are  vectors 
of  the  angular  velocity  and  acceleration  of  the  carrier,  i,  j,  and  k  are 
unit  vectors  of  system  O^XYZ,  and  p  =  xi  +  yj  +  zk;  differentiation  with 
respect  to  t  is  denoted  by  the  dot. 

Since  vectors  a^ ,  aj- ,  and  a^ ,  used  for  calculation  of  f ,  are  calculated 
in  coordinate  system  OcXYZ,  we  initially  find  its  projections 


fx 

iy 


f 

/  : 


-  r  -  (Q;,  --  Q:)  x  ^  --  x  f  2  (ClyZ  -  Qj,)}, 

-  T  -r  —  (^A-  -f  +  ilp.-z]  -f  ^  +  2  (9..x  —  q  J)], 

y  Pxy  —  -rQyjz-^  ^-x-Pi  -  z  2  (Q,.i/  — 


from  vdiich  by  formulas 

-  [jji^  --  /,,n„  ^  /-nr)  sin  7.  -r  (fx^x  -  -r  / 7.- 
/„  =  /,T,  -  f,,Xy  -  /.-r  (7.  =  arctg  f/7,  f/)) 

we  turn  to  components  fu,  fv»  snd  fw,  used  in  equations  (1). 


fa  = 
/r  =-•  - 
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Having  substituted  ordinary  derivatives  for  independent  variables  s  and 
t  in  equations  (l)-^(4)  by  partial  derivatives,  we  find  the  solving 
system  of  nonlinear  equations  of  relative  vibrations  of  the  rings • 

When  studying  the  dynamics  of  the  entire  system,  one  must  construct  the 
equations  of  motion  of  the  carrier  and  solve  them  jointly  with  the 
equations  of  motion  of  the  icings.  To  do  this,  let  us  use  the  theorems 
on  variation  of  the  momentum  Q  and  of  the  mcMnent  of  momentum  K  of  the 
carrier  in  the  forms 


di  '  ^  '  di 


^\-  -  M' 


(6) 


where  R  and  are  the  principal  vector  and  principal  moment  of 
reactions  at  the  points  of  finishing  of  flexible  rings,  caused  by  their 
elastic  deformation  upon  vibrations  of  the  system,  and  and  are  the 
principal  vector  and  principal  moment  of  all  active  forces. 

Using  equations  (l)-(6),  the  dynamics  of  the  induced  and  free  motions  of 
the  system  is  studied.  Let  us  dwell  on  the  results  of  solving  the 
problems  of  free  oscillations. 

Let  us  limit  ourselves  to  calculation  of  the  lower  values  of  the  range 
of  vibration  frequencies  and  to  construction  of  the  corresponding  forms. 
This  permits  one  to  turn  from  the  initial  continual  problem  to  the 
corresponding  discrete  problem  and  to  consider  only  a  finite  ntnnber  of 
degrees  of  the  system. 

Let  us  assume  that  the  vibrations  of  the  system  are  determined  by 
time-dependent  displacements  of  the  system  of  its  discrete  points 

where  N  is  the  number  of  discrete  points.  The 

equations  of  free  motions  of  the  system  can  then  be  represented  in 
matrix  form 


MX  KX  -  0, 


(7) 


where  X  is  the  vector  of  displacements  of  the  nodal  points  of  the 
system,  M  is  the  matrix  of  masses,  and  K  is  the  matrix  of  the  stiffness 
of  the  discrete  elements. 

Matrix  M  is  determined  by  the  mass  of  each  extended  element, 
concentrated  at  the  nodal  points. 

The  stiffness  matrix  K  can  be  determined  by  using  the  property  K  =  $’  * , 
where  $  is  a  pliancy  matrix.  The  elements  of  matrix  $  are  the 
displacements  along  the  permissible  couplings  of  the  directions  upon 
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action  of  the  corresponding  unit  concentrated  force  on  the  node  [ 6 ] . 
They  are  found  as  a  result  of  solving  the  problems  of  the  dynamic 
equilibrium  of  the  rings  upon  corresponding  steady  translational  and 
rotational  motions  of  the  system  [4]. 

Assuming  that  free  motion  of  the  system  heus  the  form  of  harmonic 
vibrations,  let  us  represent  the  solution  of  the  system  in  the  form 
X{t)  =  H  sin  at,  »diere  H  is  the  vector  of  the  form  of  vibrations  and  a 
is  the  circular  frequency  corresponding  to  it.  Substituting  this 
solution  into  (7),  we  find 


(/C  — o)-M)H  =  0. 


(8) 


The  eigenvalues  and  the  eigenvectors  Hi  are  calculated  using  the 
QR-algorithm. 


(1) 

(2) 

(3) 

(4) 

riet^eincHHe  no 

HoMep  tbopMfai 

ricpHOA,  c  (T^*) 

qacroTa,  paA/c  (<o,*) 

Vror,  noBOpora,  paA 

HanpaBJicHHio  ocu 

>  1 

1  - 

3 

1 

5 

1 

175,75 

0,0357 

— 

Ox 

2 

132,54 

0,0474 

==0,0189 

Ox 

3 

62,46 

0,1006 

— 

Ox 

4 

57.02 

0,1102 

qr^=:0,0159 

Ox 

0 

6 

40,13 

0,1566 

— 

Ov 

26,86 

0,2340 

— 

Ox 

7 

26,76 

0,2348 

(f  ^^-0,00389 

Ox 

8 

15,65 

0,4016 

cp^= — 0,0387 

Oy 

9 

14,29 

0,4398 

Ox 

10 

14,26 

0,4408 

— 

Oz 

11 

14,24 

0,4411 

fp-^ -0,000008 

Oz 

12 

14,15 

0,4439 

fp  =0,0319 

Oy 

8,96 

0,7014 

(p%0,0032 

Ox 

14 

8,95 

0.7017 

(Py=0, 00105 

Ox 

15 

6,82 

0,9207 

— 

Oz 

16 

6,37 

0,9861 

II 

O 

o 

4^* 

00 

Oz 

17 

6,333 

0,9921 

— 

Ox 

18 

6,328 

0,9930 

ir  ^,=0.0419 
.p‘^=:0,0003 

Oz 

KEY: 


1. 

Number  of  form 

4.  Angle  of  rotation,  rad 

2. 

Period,  s(Ti) 

5.  Displacement  in  direction  of 

3. 

Frequency,  rad/s  {(Ji) 

axis 

A  total  of  18  values  of  the  spectrum  of  the  natural  vibration 
frequencies  of  the  carrier-elastic  ring  system  is  calculated  by  the 
developed  method,  and  the  corresponding  forms  of  motion  are  constructed. 
The  results  of  calculation  are  presented  in  the  table.  The  ordinal 


5 


numbers  of  the  eigenvalues  are  given  in  the  first  column,  the  length  of 
the  vibration  period  Ti  is  given  in  the  second  column,  the  corresponding 
values  of  the  circular  frequencies  Ui  are  given  in  the  third  column,  the 
amplitude  values  of  the  angles  of  rotation  of  the  carrier  with  respect 
to  the  corresponding  axis  are  given  in  the  fourth  column,  and  the  axis, 
in  the  direction  of  which  the  greatest  displacement  of  the  typical  point 
of  the  ring  is  realized,  is  indicated  in  the  fifth  column.  The  forms  of 
vibrations  of  the  rings  are  normalized  so  that  the  amplitude  value  of 
the  maximum  displacement  is  equal  to  one. 

The  axonometric  images  of  the  first  six  forms  of  natural  vibrations  are 
shown  in  Figure  1  with  respect  to  a  fixed  coordinate  system.  Two 
coordinate  systems — a  fixed  system  OXYZ  and  coordinate  system  Qxyz, 
rigidly  bound  to  the  carrier,  rotated  during  its  vibrations  with  respect 
to  the  fixed  system — are  also  presented  in  it.  The  movable  coordinate 
system  does  not  happen  in  cases  when  the  rotational  motion  of  the  body 
is  absent.  The  forms  of  vibrations  of  the  rings  are  normalized  in  all 
figures  so  that  the  value  of  the  maximum  relative  displacement  is  equal 
to  1/5  of  the  diameter  of  the  ring. 

As  calculations  showed,  the  vibration  modes  that  permit  rotation  of  the 
carrier  have  higher  frequencies,  compared  to  the  descending  modes  of 
natural  vibrations  of  the  rings,  rigidly  fastened  to  a  fixed  support. 

Upon  study  of  the  dynamic  and  strength  characteristics  and  also  of  the 
controllability  of  the  cargo  spacecraft  with  unfolded  large  annular 
structures,  the  following  problems  were  solved: 

1 )  experimental  determination  of  the  frequency  and  decrement  of  the  main 
tone  of  vibrations  of  the  structures; 

2)  a  check  of  the  capability  of  the  structures  of  restoring  a  given  form 
after  the  cargo  spacecraft  completed  turns  and  correction  pulses; 

3)  a  check  of  the  results  of  preliminary  calculations  of  the  structural 
djTiamics  and  applicability  of  the  calculating  models  to  be  used. 

The  most  complete  information,  including  a  television  image  of  the 
vibrations  of  one  of  the  annular  structures,  transmitted  from  the  cargo 
spacecraft  (Figure  2),  during  the  experiment  was  found  in  two  sessions, 
in  viiich  the  cargo  spacecraft  accordingly  performed  a  correcting  pulse 
Av  =  0.6  m/s  and  a  programmed  yawing  turn  by  angle  jp  =  120  at  angular 
velocity  u  =  0.67”/s. 

Synmetrical  vibrations  of  the  structures  were  excited  with  period  190  s 
during  the  correcting  pulse.  The  points  of  the  rings  most  remote  from 
the  carrier  were  deflected  from  the  neutral  position  by  11  m.  piese 
values  of  the  period  and  amplitude  of  vibrations  essentially  coincided 
with  the  results  of  earlier  calculations.  It  should  be  noted  that  the 
actual  period  of  vibrations  was  determined  immediately  from  three 
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sources:  by  a  television  image  frcm  the  cargo  spacecraft,  by  a  report 
of  the  crew  of  the  "Mir"  spacecraft  and  by  telemetry  data. 


Figure  1 

The  dissipative  properties  of  the  large  annular  structures  were  studied 
upon  analysis  of  the  vibrations  caused  by  the  correcting  pulse,  which 
are  difficult  to  determine  by  calculation  or  by  ground  tests.  Observing 
these  vibrations  directly,  the  crew  of  the  "Mir"  orbital  station  noted 
slow  damping  of  them,  while  preliminary  analysis  of  the  telemetry  data 
permitted  an  estimation  of  the  logarithmic  decrement  of  S  -  0.05. 


Figure  3 


When  considering  the  methods  of  improving  the  controllability  of  a 
spacecraft  with  large  unfolded  annular  structures  and  taking  into 
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account  the  weak  danipinS  of  vibration  of  these  structures ,  main 
attention  should  be  devoted  to  dynamic  suppression  of  vibrations, 
excited  during  maneuvers  of  the  spacecraft,  rather  than  to  passive 
damping.  Thus,  the  dynamic  effect  of  essentially  total  instant^eous 
correction  of  vibrations  of  the  structure  as  a  result  of  selecting  the 
moment  of  stopping  the  rotation  of  the  spacecraft  can  be  specifically 
used  as  one  of  the  principles  of  selecting  its  optimal  control. 


Figure  4 

The  dynamic  effect  of  essentially  total  kinematic  suppression  of 
vibrations  was  demonstrated  vdien  the  Progress-40  spacecraft  performed  a 
programmed  turn,  and  the  angle  of  rotation  of  (?  =  120  was  previously 
calculated  on  the  basis  of  the  measured  value  of  the  actual  period  of 
the  main  tone  of  vibrations  of  the  annular  structures— T  =  190  s^and  of 
the  ncaninal  rotational  velocity  of  the  cargo  spacecraft  u  =  0.67  /s, 
with  the  calculation  that  the  structure  performed  a  single  vibration 
during  the  turn  of  the  cargo  spacecraft. 

It  was  obvious  from  the  television  image  from  the  Progress-40  spacecraft 
(Figures  2-4)  that  the  annular  structure  performed  one  vibration  during 
the  turn,  and  that  it  essentially  stopped  in  the  neutral  position  at  the 
moment  of  completion  of  the  turn  and  the  beginning  of  stabilization  of 
the  cargo  spacecraft.  The  theoretical  conclusion  of  the  possibility  of 
rapid  suppression  of  vibrations,  having  the  value  of  the  optimal  control 
of  the  dynamics  of  large  structures  vdien  working  out  the  methods,  was 
thus  confirmed. 
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Comparison  of  the  results  of  calculation  and  data  of  full-scale  tests 
indicates  that  they  are  in  satisfactory  agreanent. 
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ON  ERRORS  OF  SOLID-STATE  WAVE  GYROSCOPES  DURING  PROGRESSIVE  VIBRATION  OF 
BASE 

917F0159A  Kiev  DOKLADY  AKADEMII  NAUK  UKRAINSKOY  SSR:  SERIYA  A  — 
FIZIKO-MATEMATICHESKIYE  I  TEKHNICHESKIYE  NAUKI  in  Russian  No  11,  1990 
(manuscript  received  4  Jun  90)  pp  25-28 
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[Text]  The  effect  of  progressive  and  angular  vibration  on  the  motion  of 
traditional  mechanical  gyroscopes  has  now  been  studied.  However,  the 
problem  of  the  effect  of  vibration  of  the  base  on  the  accuracy  of  the 
new  inertial  infonnation  meter— the  solid-state  wave  gyroscope  (TVG) , 
vdiich  is  of  considerable  interest,  has  been  insufficiently  studied 
[1-4]. 

The  nonlinear  dynamics  of  a  solid-state  wave  gyroscope  with  imperfect 
cavity  on  a  progressively  vibrating  base  is  considered  in  this  paper. 

The  state  of  the  object  is  given  by  generalized  coordinates  pn(t)  and 
qn(t) ,  corresponding  to  conjugate  natural  forms  of  vibrations  of  a 

hemisfiierical  cavity  l/hn(0)cos /i<p,  t/ftn(0) sin ncp  (k=l,  2,  3),,  which  are 

introduced  by  the  approximation 


Wi  =r  VUi,,  (p„  COS  «fp  -f  q^.  sin  nq), 

n 

^2  =  V  1^271  (Pn  sin  AZff-  ^  t/,,  C()c>  /n; 

n 

COS  «rp  --  q^  sin  n(f>). 


where  angle  p  is  counted  in  the  circumferential  direction  and  angle  0  is 
formed  by  the  axis  of  symmetry  and  by  the  internal  normal.  Application 
of  the  Bubnov-Galerkin  procedure  to  the  nonlinear  equations  of  motion  of 
the  cavity  on  a  vibrating  base  [4]  with  regard  to  representation  of 
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defects  of  a  hemispiierical  dome  in  the  form  of  expansions  into  a  Fourier 
series  by  circumferential  angle  tp  [5]  results  in  ordinary  differential 
equations  in  generalized  coordinates  that  describe  nonlinear  vibrations 
of  an  imperfect  cavity  on  a  progressively  vibrating  base.  Let  us 
consider  resonant  vibrations  idien  the  external  perturbation  frequency  x 
is  a  multiple  of  the  working  frequency  ^2  (the  case  of  excitation  of  a 
second  pair  of  natural  forms  of  vibrations  in  the  cavity,  n  =  2,  is  of 
interest  for  applied  problems  [1,  6]),  assuming  that  there  is  no 
internal  resonance  in  the  system.  Equations  that  describe  the  dynamics 
of  the  main  form  have  the  form 


(mo  +  e/nii)  q,  +  +  (^o  +  P2  +  C12P2  + 

+  e  (//oPo?2  -r  +  0 

(mo  -f  e/r.22)  'p2  +  e«2i<72  —  2eA'o^272  -r  e/-oP2  +  (^o  +  ^<^22)  -t-  C21P2  t 

4-  e  (H0P0P2  +  -r  e/2  +  0  (e-)  =  0. 


(1) 


where 


=  rriiWox  +  ov  /'•:  —  .  /Ho^ov  +  /tt.U 


Here  moi  co  t  Ko  and  .lo  are  the  coefficients  of  inertia  and  stiffness, 
and  gyroscopic  and  damping  coefficients  of  an  ideal  cavity, 
respectively,  ft  is  the  angular  rotational  velocity  of  the  cavity  about 
its  own  axis  of  symmetry  (it  is  considered  a  slowly  variable  time 

function,  i.e.,  we  disregard  terms  with  ft  and  ft^),  m^j  and  Cij  (i,  j  = 

1,  2)  are  perturbations  of  the  coefficients  of  inertia  and  stiffness  due 
to  errors  in  manufacture  of  the  hemispherical  cavity  (the  fourth 
harmonics  of  expansion  of  the 

defects  into  a  Fourier  series  [5]),  and  =  012  —  021,  rrin— — ^22, 

cii=  — C22:  ^7^’.  (^=1.  2,  3)  are  perturbations  of  the  inertial 

characteristics  of  the  cavity,  m3  and  ma^  correspond  to  the  second 
ccffnponent  of  the  defect  in  Fourier  expansions,  mi,  mi*,  m2,  and  m2*  are 

the  first  and  third  harmonic,  =  — //Aocosv/,  Wqy— — 

Woz=  — x^Zo  cos  xt  are  projections  of  the  vector  of  linear  accelerations 

of  the  base  onto  axes  X,  Y,  and  Z  of  a  coordinate  system  rigidly  bound 
to  pole  O  of  the  cavity,  and  axis  OZ  coincides  with  the  axis  of  symmetry 
of  the  cavity,  qi ,  pi ,  and  po  are  generalized  coordinates  that 
characterize  the  motion  of  a  hemispherical  cavity  as  a  solid 


cos  yi,  p,  =  g\Y,  cos  yt.  Po  =  yt,  g,  =  const, 

=  const,  =  const, 
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where  Xo,  Yo»  Zo  and  X  are  the  amplitudes  and  frequency  of  external 
p)erturbation ,  Ho»  Hit  H21  and  H3  are  constant  coefficients  with 
nonlinear  terms  1  and  e  is  a  small  parameter. 

Let  us  consider  two  important  cases  of  vibration  of  the  base. 

Longitudinal  vibration.  The  base  moves  along  the  axis  of  symmetry  of 
the  cavity  (Z  =  Zo  cos  Analysis  of  equations  (1)  shows  that  the 

following  resonance  situations  are  possible  in  the  system:  a)  ^  =  ^2  + 
+  €(r,  idiere  (r  is  a  detuning  parameter.  In  this  case,  longitudinal 
vibration  of  the  base  results  in  excitation  of  "spurious"  induced 
vibrations  in  the  cavity.  In  zero  approximation  of  the  perturbation 
method  [7],  solution  of  the  system  of  equations  (1)  has  the  form 


qW  =  Q  (Tj)  exp  {'mJo)  +  Q  (Ti)  exp  (— 
p(0)  ^  p  (7  j  exp  (( co„r„)  +  P{Ti)  exp  (—  tco.T „), 


where  7'„  =  8’‘/;  Q,  Q  and  P,  P  are  complex  conjugate  amplitudes.  Let  us 
represent  the  complex  amplitudes  in  exponential  form 


Q  =  a  (Ti)  exp  {t  [a  (T 1)  +  aT J}, 
P  =  hiT,)exp[ilf,(T,)-^aT,i}. 


Substitution  of  variables  by  formulas  b  +  ia  —  C exp (2i'F),  p  —  <2=4),  a  +  p  = 

=  7  permits  us  to  write  the  condition  of  the  absence  of  secular  terms 
in  the  solutions  of  the  first  approximation  q^  •  *  2  p^**2  i*^  new 

variables  C,  f,  p,  and  7.  A  real  instrument  has  a  control  system  that 
maintains  vibrations  of  the  cavity  in  the  form  of  a  standing  wave, 
since,  as  shown  in  [8],  it  assumes  fulfillment  of  the  following 

conditions:  ip  ^  Taking  the  foregoing  into  account, 

we  find  a  formula  that  determines  the  precession  rate  of  the  standing 
wave  during  functioning  of  a  solid-state  gyroscope  due  to  the  effect  of 
longitudinal  vibration 


T'  = 


K,P.  , 

2m, 


%'Z,  sin  (y/2) 


{m,  cos  2'F  —  m’  sin  2¥). 


(2) 
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Here  (•••)=  *^1^,---;  C  =  (6- -f- 'F  are  the  vibration  amplitude  and  angle 

ul  I 

of  orientation  of  the  wave  pattern,  respectively,  and  =  is 

the  natural  frequency  of  an  ideal  cavity.  It  follows  frcsn  formula  (2) 
that  the  integrating  effect  of  a  solid-state  wave  gyroscope  will  occur 
only  at  angular  velocities  ft,  which  exceed  the  threshold  value 

^min  =  AmZo(o.,/2A'„C.  Here  Am  =  (m|  -j-  %  =  ®2- 

b)  X  -  24^2  +  2e(r.  Longitudinal  vibration  excites  "spurious"  parametric 
vibrations.  The  equation  that  determines  angle  f  has  the  form 

T  =  -K^/2m,,  (3) 


It  follows  from  equation  (3)  that  longitudinal  vibration  in  the  given 
case  does  not  result  in  an  additional  "spurious"  component  in  the  wave 
pattern. 

Transverse  vibration.  The  base  moves  in  a  plane  perpendicular  to  the 
axis  of  symmetry  of  the  cavity  (A=AoCos"/^,  y=Focosx/).  The  following 
resonant  situations  are  possible  in  the  system: 


a)  %  =  -r  ec7. 


1 

b)  X  ==  — 


2 


ea. 


In  both  cases,  transverse  vibration  of  the  base  results  in  excitation  of 
"spurious"  induced  vibrations  of  the  cavity.  Expressions  that  determine 
the  position  of  the  wave  pattern  as  a  function  of  the  anguleu"  rotational 
velocity  and  vibration  parameters  of  the  base  have  the  form 


a) 


b) 


'F  —  (2m„)  '  A(,Q  +  (4m„(u„Q  '  m„Vg)cos2T — 

—  (m*A"j,  -f  m'Kj)  sin  2T]  sin  {yl2), 

\Hj^g^g\X,,\\c.os,2W — 

-  ~  sin  2^^  sin  {y!2). 


In  the  given  case,  we  will  also  have  some  threshold  value  of  the  angular 
rotational  velocity  of  the  base,  below  which  there  is  no  integrating 
effect. 

It  follows  from  the  studies  that  there  are  several  mechanisms  of  errors 
in  the  operation  of  the  instrument  upon  vibration  of  the  base: 
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excitation  of  additional  "spurious"  vibrations  of  the  cavity  that  leewi 
to  distortion  of  the  readings  of  the  instrument  due  to  the  presence  of 
errors  in  manufacture  of  the  cavity; 

excitation  of  additional  "spurious"  vibrations  due  to  the  presence  of 
nonlinear  coupling  (between  generalized  coordinates)  that  characterize 
intermodal  interaction. 

The  results  of  the  studies  can  be  used  in  determining  the  dimensions  of 
the  zones  of  absence  of  the  integrating  effect  of  a  solid-state  wave 
gyroscope,  the  nominal  value  of  the  vibration  amplitude  of  the  cavity, 
and  also  when  working  out  the  requirements  on  the  precision  of 
manufacture  of  the  cavity  and  on  the  vibration  protection  system. 
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[Article  by  V.  A.  Zebergs,  L.  L.  Eihmanis,  and  I.  Ya.  Niedrite] 

[Text]  The  problems  of  increasing  the  efficiency  of  utilizing  fuel  and 
energy  resources  are  acquiring  ever  more  timely  value  in  all  sectors  of 
the  national  economy.  Problems  of  energy  conservation  will  determine  to 
a  great  extent  in  the  future  the  development  of  industry,  transport, 
agriculture  and  the  communal-everyday  sector.  The  problem  has  special 
significance  for  Latvia  ijnder  conditions  of  transition  to  independent 
forms  of  management  and  republic  cost-accounting,  since  85  percent  of 
the  fuel  and  50  percent  of  the  electric  power  are  delivered  from  other 
regions.  Along  with  this,  the  problem  of  adequate  and  real  nature- 
conserving  activity  is  valid.  Interaction  of  facilities  of  the  fuel  and 
energy  complex  with  the  environment  is  an  exceptionally  important 
problem.  On  the  one  hand,  the  value  of  an  uninterrupted  effective 
energy  supply  of  all  sectors  of  the  national  economy  is  very  high,  and 
on  the  other  hand,  the  fuel  and  energy  complex  has  the  most  significant 
effect  on  the  human  habitat.  Only  centralized  power  enterprises 
discharge  more  than  40,000  tons  of  noxious  materials  annually  in  the 
atmosphere  in  Latvia,  which  comprises  21  percent  of  the  total  discharge 
in  the  republic. 

The  increasing  pollution  of  the  atmosphere  is  assuming  a  global  nature 
today.  Dischairges  of  noxious  and  harmful  pollutants  into  the  atmosphere 
by  many  sectors  of  industry  have  a  deleterious  effect  on  living 
orgeinisms,  destroy  the  ozone  layer,  and  contribute  to  the  occurrence  of 
acid  rain.  The  reality  of  the  greenhouse  effect  occurring  is  a 
consequence  of  variation  of  the  structure  of  radiaint  heat  transfer  as  a 
result  of  intensive  utilization  of  thermal  power  systems.  These 
phenomena  may  affect  the  entire  country  and  even  vdiole  continents. 
Although  scientists  have  different  opinions  on  the  time  intervals  and 
intensity  of  the  ecological  changes,  nevertheless  they  are  still 
significant,  which  requires  urgent  study. 

A  session  of  the  General  Assembly  of  the  Academy  of  Sciences,  devoted  to 
discussion  of  the  project  of  the  main  propositions  of  the  Program  of 
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Biospheric  and  Ecological  Research  of  the  USSR  Academy  of  Sciences  for 
the  period  up  to  2015  in  December  1988  noted  that  we  are  in  the  phase  of 
anthropogenic  effects  on  nature,  when  the  processes  of  ecological 
degradation  are  even  more  local  in  nature,  but  global  trends,  hazardous 
for  the  biosphere  as  a  whole,  are  already  being  detected  on  this 
background.  This  is  the  moment  when  the  process  of  irreversible 
ecological  changes  can  be  stopped  now  through  the  efforts  of  people. 

Outside  of  health  and  well-being,  all  our  social  goals  and  political  and 
ecological  ideals  are  impractical.  Alternative  versions  of  solving 
socioecological  and  technical  problems  must  be  worked  out.  Academician 
P.  L.  Kapitsa  convincingly  demonstrated  the  interdependence  of  the 
growth  of  the  gross  national  product  and  energy  production  for  all 
countries  of  the  world.  This  means  that  if  there  is  an  increase  of 
energy,  there  will  also  be  an  increase  of  the  national  product. 

The  time  has  come  to  combine  the  means  of  survival  into  a  single 
complex,  optimal  in  the  given  ecological  situation  [11.  This  is  a  most 
important  scientific  task  both  on  a  global  scale  and  for  each  country. 

The  Physics  and  Power  Institute  of  the  Latvian  Academy  of  Sciences  is 
participating  in  developnent  of  the  section  '’Power  engineering"  of  the 
general  academic  program  of  biospheric  and  ecological  research  for  the 
period  up  to  2015.  The  topic  being  studied  is  a  method  of  estimating 
the  loss  due  to  the  effects  of  energy  facilities.  Problems  of 
estimating  the  ecological  impact  of  energy-conserving  measures,  vAiich 
include  analysis  of  the  engineering  solutions  of  energy  conservation  in 
industry  and  ecological  zoning  and  dipersion  of  energy  resources, 
structural  shifts  in  industry,  and  their  energy  and  ecological 
consequences,  are  being  considered  in  the  first  phase  of  these  studies. 
The  studies  are  directed  toward  discovering  the  characteristic  features 
of  energy  conservation  with  different  layouts  of  energy  conservation  of 
enterprises,  a  rational  structure  of  industry  and  with  specific 
principles  of  estimating  the  ecological  impact  of  energy-conserving 
plants,  now  used  in  Latvia. 

Determination  of  the  principles  of  ecological  estimation  of  the 
engineering  solutions  of  energy  conservation  in  our  republic  is  based  on 
analysis  of  the  characteristic  features  of  natural  conditions, 
functioning  of  the  natural  and  technical  environment,  and  of  the  state 
of  the  ecoosystems  of  the  given  territory.  Thus,  optimization 
calculations  are  made  with  regard  to  ecological  factors  under  specific 
conditions  of  a  corresponding  region,  but  the  principal  approach  to  the 
problem  is  identical  for  all  economic  regions. 

Analysis  of  the  functioning  of  the  territorial  economic-ecological 
system  of  the  Latvian  Republic  shows  that  its  intensity  determines  the 
consequences  of  the  ecological  concept,  implemented  during  the  previous 
5  years,  which  was  constructed  on  the  basis  of  hypertrophic  development 
of  industry,  oriented  toward  extensive  growth  factors,  and  was  not 
accompanied  by  a  search  for  effective  ways  of  reducing  the  energy  and 
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resource  capacity,  and  solution  of  environment  protection  problems.  The 
main  industrial  centers  at  present — cities  and  urban  agglomerations — 
have  become  areas  of  extensive  negative  transformation  of  the 
environment,  the  effect  of  which  is  very  significant  not  only  on  the 
nature  of  adjacent  territories,  but  also  on  the  vast  regions  of  river 
basins,  and  coastal  waters  of  the  sea.  One  must  sometimes  talk  about 
the  symptoms  of  a  crisis  state,  when  the  parameters  of  the  environment 
are  approaching  the  permissible  limits  of  changes,  and  suppression  of 
reproductive  capacities  and  exhaustion  of  the  ecological  capacity  of  the 
territory  have  been  noted;  therefore,  intensification  of  the 
environmental  protection  trend  now  conducted  outside  economic  reform  is 
necessary,  the  ecology  aspect  of  cost  accounting  is  required,  and  this 
must  be  transformed  to  a  reliable  barrier  on  the  path  toward  wasteful 
utilization  of  nature. 

The  current  phase  of  development  of  the  economy  is  characterized  by  a 
stable  increase  of  the  expenditures  of  social  labor,  embedded  in 
environmental  protection  and  in  reproduction  of  natural  resources.  A 
growth  of  the  scales  of  ecological  investment  increases  the 
responsibility  of  each  planning  solution  in  this  field,  and  imposes 
higher  requirements  on  their  effectiveness.  At  the  same  time,  society 
needs  optimal  development  of  TEK  [not  further  identified]  as  a  necessary 
condition  of  enhancing  the  life  level  and  national  income.  An  objective 
need  arises  in  this  regard  to  establish  interrelationships  between  the 
results  of  economic  activity  and  its  ecological  consequences. 
Substantiation  of  the  economy  of  production  is  an  inseparable  part  of 
the  management  system  that  affects  the  selection  of  priorities  in 
support  of  the  national  economy  udth  natural  and  energy  resources.  The 
main  causes  of  a  negative  effect  of  the  national  economy  on  the 
environment  are  imperfection  of  technology,  an  increase  of  the  scales  of 
production,  and  its  concentration  on  a  limited  territory.  The  increase 
of  energy  consumption  can  be  covered  by  30-40  percent  in  the  Latvian 
republic  by  introduction  of  energy-conservation  technologies  and 
scientifically  substantiated  standards  of  energy  consumption  [2]. 

Life  persistently  requires  extensive  analysis  of  the  energy  consumption 
of  manufactured  products.  High  losses  of  energy  resources  are  formed 
due  to  the  fact  that  we  do  much  that  is  unnecessary,  and  even  in 
weighted  form.  The  unsubstantiated  high  materials  consumption  is 
accompanied  by  impermissibly  high  energy  consumption.  Recovery  of 
secondary  resources  that  yields  significant  energy  conservation  has  not 
become  widespread.  The  high  degree  of  depletion  of  production  funds  in 
industry,  which  comprises  more  than  47  percent,  also  results  in  rational 
consumption  of  energy  in  our  republic.  This  part  of  the  morally  and 
physically  obsolete  production  funds  is  also  a  significant  source  of 
environmental  pollution. 

The  fact  that  tenfold  more  funds  are  expended  in  the  country  on  routine 
maintenance  of  operated  equipment  than  is  required  for  replacement  of  it 
indicates  the  redundancy  of  this  practice  [3]. 
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utilization  of  secondary  energy  resources  can  reach  more  than  46  percent 
in  the  future  structure  of  measures  on  conservation  of  energy  resources, 
while  the  use  of  energy— conser\'^ing  technology  in  industry  and  in 
agriculture  may  reach  approximately  38  jjercent  [21. 

An  experimental  lot  of  facade  heat  regulators  is  functioning 
successfully  in  Riga.  The  instruments,  developed  at  FEI,  contribute  to 
conservation  of  energy  resources  and  prevent  overconsumption  of  thennal 
energy.  The  heat  loss  factor  of  new  multistory  buildings  is  0.42-0.47 
kilocalories  per  square  meter,  which  is  25-30  percent  higher  than  the 
optimal  value. 

Taking  into  account  the  ecological  impact  is  of  considerable 
significance  in  economic  estimation  of  energy— conserving  measures  and  in 
determination  of  the  priority  of  implementing  them  with  limited  monetary 
and  material  funds.  Analysis  of  the  ecological  impact  is  related  to 
estimation  of  the  consequences  of  discharges  of  energy  sources  that 
supply  energy-conserving  plants  compared  to  treiditional  plants,  which 
utilize  a  closed  form  of  fuel  (in  the  presence  of  scrubbing  structures 
that  provide  the  PDV  [maximum  permissible  discharge]  at  the  required  PDK 
[maximum  permissible  concentration] ) .  The  consumption  of  energy 
resources  on  the  energy  supply  of  the  drive  of  energy-conserving  plants 
for  use  of  secondary  thermal  resources  is  insignificant,  and  sometimes 
(utilization  of  TNU  heat  pumps)  may  reach  30  percent  or  more  of  the 
extracted  alternative  energy  resources.  Therefore,  the  ecological 
impact  is  determined  by  the  difference  in  the  amount  of  fuel  burned  for 
comparable  variants  (the  drive  of  plants  that  conserve  energy  or 
generate  the  same  amount  of  energy  by  the  traditional  method)  and  by  the 
loss  due  to  discharges.  Estimation  of  the  loss  is  determination  of  the 
expenditures,  related  to  usefully  generated  (conserved)  energy  at 
comparable  energy  plants  ( per  ton  of  conventional  fuel ) .  The  economic 
impact  of  energy-conserving  measures  can  generally  be  formulated  with 
regard  to  the  ecological  loss  due  to  discharges  as 
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where  ^3'  .  is  the  specific  economic  impact  of  energy-conserving 

vC 

measures  with  regard  to  the  ecological  impact  (rubles/t  of  conventional 
fuel).  As  is  the  difference  of  the  specific  expenditures,  dependent  on 

capital  investments  on  a  traditional  power  plant  operating  on  the  close 
fuel  cycle  and  on  an  energy-conserving  plant,  A3^  is  the  difference  of 

the  specific  expenditures,  dependent  on  operating  expenditures,  A3,j.  is 

the  difference  of  the  specific  expenditures  for  fuel  and  energy  for 
comparable  variants ,  and  Ay  is  the  difference  in  the  loss  due  to 
discharges . 
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The  difference  of  specific  fuel  and  energy  expenditures  characterizes 
the  conservation  of  energy  resources  and  can  be  determined  as 
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which  occurs  upon  compilation  of  alternative  variants, 
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where  B  and  B  are  fuel  consumption  for  a  traditional  energy  plant 
3  T  3C 

that  utilizes  a  complete  type  of  fuel  and  an  energy-conserving  plant 
(tons  of  conventional  fuel),  and  are  the  useful  energy 

consumption  (tons  of  complete  conventional  fuel),  and  are  the 

energy  utilization  factors  from  0.8  to  0.5  for  traditional  heat 

supply  plants,  while  ^  =  1.0  to  2.0  or  more  for  energy-conserving 

plants,  for  example,  for  TNU,  due  to  the  use  of  ambient  thermal  energy), 
and  z^^  are  the  specific  expenditures  for  a  complete  type  of  fuel 

( rubles/ton  of  conventional  fuel ) . 

The  difference  in  damage  due  to  discharges  Ay  (upon  observation  of  PDV 
and  PDK)  is  dependent  on  the  difference  in  the  amount  of  burned  fuel  per 
unit  of  useful  energy  and  of  the  specific  damage  due  to  discharges  in 
the  compared  variants ,  i . e . , 
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where  Ay  is  the  difference  in  damage  due  to  discharges  (rubles/t  of 
complete  conventional  fuel)  and  y^  is  the  specific  damage  due  to 

discharges  per  unit  of  burned  complete  type  of  fuel  (rubles/t  of 
conventional  fuel ) . 

When  using  different  types  of  fuel  to  drive  an  energy-conserving  plant 
and  in  a  traditional  plant.  Ay  is  determined  as 
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where  yr,  and  y_  are  the  specific  damage  due  to  discharges  from 

combustion  of  complete  type  of  fuel  and  of  fuel  utilized  to  generate 
energy  to  drive  an  energy-conserving  plant  (rubles/t  of  conventional 
fuel),  respectively. 

Selecting  the  type  of  complete  fuel  and  energy  to  drive  energy- 
conserving  plants  is  related  to  optimization  of  the  fuel-energy  balance 
of  the  region,  the  structure  of  which  can  be  changed  significantly  upon 
implementation  of  wide-scale  measures  on  conservation  of  energy 
resources  (gassification  of  a  region).  To  determine  the  volume  of 
conservation  of  energy  resources  and  of  the  amount  of  complete  type  of 
fuel,  one  must  also  resort  to  using  optimization  models. 

A  method  that  determines  the  damage  due  to  discharges  Ay  is  rather 
complex  and  requires  further  studies.  However,  this  is  a  separate 
direction  of  the  studies,  related  to  participation  of  scientists  of 
different  specialties  on  problems  of  ecology. 

Taking  into  account  the  additional  economic  effect  of  energy-conserving 
measures  in  a  region  significantly  affects  both  their  ranking  according 
to  ecological  impact,  and  the  volume  of  feasible  implementation  of 
energy-conserving  measures  for  the  futtire. 
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[Text]  General  Characteristic 

The  program  system  (PK)  is  designed  to  study  the  main  directions  and 
principles  of  development  of  energy  associations  (EO).  The  most  typical 
problems  to  be  solved  by  using  the  ORS  E  program  system  are:  1) 
selecting  the  variants  of  development  of  systems-forming  networks  of 
energy  associations  at  existing  voltage  by  introducing  new  electric 
power  transmission  lines  or  by  renovation  of  them;  2)  substantiation  of 
the  introduction  of  electric  power  transmission  at  the  highest  voltage 
with  regard  to  different  times  of  implementation  of  the  corresponding 
measures  (for  example,  with  regard  to  temporary  operation  of  new 
electric  power  transmission  at  reduced  voltage ) :  3 )  selection  of  the 
directions  and  v'oltage  of  new  electric  {xjwer  transmission,  designed  to 
generate  the  power  of  new  power  plants ;  4 )  refining  the  location  of  new 
power  plants  and  substantiation  of  the  plan  of  introduction  of  the 
generating  capacities  with  regard  to  expenditures  for  development  of  the 
basic  network  of  the  electric  power  system;  5)  substantiating  the 
feasibility  of  introducing  DC  electric  power  transmission. 

It  was  necessary  to  solve  a  number  of  methodical  problems  during 
development  of  the  program  system.  They  were  related  to  the  need  to 
increase  the  number  of  variables  and  of  the  number  of  calculating  modes, 
to  optimization  of  the  flow  distribution  in  a  complex  network,  which 
combines  DC  and  AC  electric  transmission,  with  regard  to  limitations  on 
capital  investments.  Using  the  ORS  E  progratm  system,  one  can  simulate 
the  process  of  development  of  energy  associations  in  different 
situations.  Alternative  measures  on  development  of  the  network  can  be 
selected  by  using  the  program  system,  their  efficiency  for  the  entire 
range  of  anticipated  future  conditions  can  be  determined  with  regard  to 
space  and  time  relationships,  one  can  determine  the  factors  that  tiave  a 
decisive  influence  on  the  optimal  process  of  development ,  and  one  can 
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determine  which  of  the  alternative  measures  and  Vvhen  they  should  be 
implemented  to  guarantee  aui  optimal  \^lue  of  the  optimization  criterion. 
Integral  criteria  for  a  sufficiently  long  calculating  period  are  used 
viien  compiling  the  strategic  variants  of  development.  This  is  related 
to  the  fact  that  the  relative  effectiveness  for  the  more  promising 
strategies  during  the  initial  period  is  lower .  compared  to  extensive 
variants  of  development.  The  number  of  compared  versions  of  development 
is  so  great  that  substantiation  of  the  solution  is  impossible  without 
special  program  systems  for  dynamic  optimization  of  the  development  of 
electric  networks.  Due  to  the  uncertainty  of  information, 
substantiation  and  decision-making  are  carried  out  on  a  dynamic  plan; 
as  the  initial  information  is  refined,  the  solution  is  also  refined. 

The  conditions  of  development  includes  information  about  the  existing 
electric  network,  about  the  generating  capacities  and  energy  consumers, 
about  future  loads,  about  the  calculating  pjeriod,  about  alternative 
future  measures,  about  restrictions,  and  about  future  technical  and 
economic  indicators. 

The  concept  of  development  step  is  used  to  describe  the  calculating 
p)eriod  in  optimization  problems.  This  is  related  to  the  fact  that 
continuous  progress  of  development  of  a  system  over  time  in  practical 
dynamic  calculation  is  always  replaced  by  a  discrete  process,  "pie  time 
segment,  during  which  the  state  of  the  system  does  not  change,  is 
considered  a  development  step.  The  ORS  E  program  system  solves  problems 
with  number  of  steps  up  to  10. 

The  inpjut  information  for  the  ORS  E  dynamic  program  system  is 
information  about  the  calculating  plan  of  the  network  and  the  calculated 
op)erating  conditions  of  the  network  during  the  calculating  pjeriod,  and 
also  about  the  set  of  possible  measures  on  development  of  the  network, 
representing  the  area  of  optimization  itself .  The  sp>atial  relations  and 
area  of  search  for  the  optimal  variant  according  to  development  steps 
are  simulated  on  the  basis  of  the  calculating  plan  of  the  network 
object.  The  main  concepts  that  characterize  the  calculating  plan  are  as 
follows : 

1)  Notes.  The  combination  of  nodes  of  the  network  facility  to  be 
optimized  is  divided  into  groups  ( subsets ) ,  which  are  used  to 
differentiate  the  technical  and  economic  indicators.  The  nodes  of  this 
group  are  characterized  by  identical  technical  and  economic  indicators 
to  compjensate  for  the  energy  losses  in  the  elements  of  energy 
associations .  It  is  assumed  in  the  ORS  E  when  calculating  the  technical 
and  econcxnic  indica.tors  that  energy  losses  in  the  elements  of  energj^ 
associations  are  covered  in  the  outer  nodes  of  these  elements,  i.e.,  in 
place. 

2)  Elements  of  network.  Consideration  of  the  following  typjes  of 
elements  of  the  network  is  p)ermissible;  a)  of  an  AC  electric 
transmission  line;  b)  of  a  DC  electric  transmission  line;  c)  of  a  power 
plant  (generators).  The  main  notation  of  the  element  in  the  primary 
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calculating  scheme  are  the  numbers  of  its  outer  nodes.  Transformers  are 
not  directly  shown  in  the  calculating  scheme.  Lines  with  different 
rates  of  voltages  may  arrive  at  one  node.  The  main  technical  and 
economic  characteristic  of  the  element  of  an  energy  association  is  the 
number  of  the  set  of  technical  and  economic  data.  Standardization  of 
elements  is  introduced  by  using  the  indicated  number.  However, 
additional  economic  indicators  (constant  components  of  capital 
investments)  can  be  used,  by  which  individual  features  of  separate 
elements  are  taken  into  account.  These  indicators  are  also  used  to  take 
into  account  the  economic  indicators  of  transformer  and  converter 
substations . 

3)  Branches.  The  elements  of  the  network  are  reflected  in  the 
calculating  plan  by  branches.  A  branch  is  characterized  by  data  of  a 
specific  element  (AC  or  DC  electric  power  transmission  line,  power 
plants)  and  by  the  number  of  parallel  elements.  Thus,  a  branch  may  also 
reflect  several  parallel  LEP  or  generators  with  identical  parameters. 

The  primary  scheme  should  contain  both  existing  and  future  elements  of 
the  network.  Thus,  branches  can  be  divided  into  existing  and  future 
branches . 

The  number  of  elements  according  to  developnent  steps  t  can  be 
increased.  Parallel  branches  that  reflect  the  elements  of  the  network 
with  different  characteristics  are  permitted  in  the  scheme. 

The  connectedness  of  the  calculating  scheme  is  checked  automatically  by 
the  programs  of  the  system.  The  layout  of  the  network  is  assumed 
connected  if  all  its  nodes,  which  can  have  a  load,  are  connected  to  each 
other  by  existing  or  future  branches.  The  characteristics  which  occur 
upon  display  of  the  electric  network  for  dynamic  optimization 
calculations  include,  for  example,  prohibition  of  elements  with  zero 
resistance  (switches).  One  can  require  that  one  object  be  represented 
by  several  branches  to  display  these  measures.  The  calculating  schemes 
of  electric  networks,  compiled  for  dynamic  optimization  calculations, 
may  have  100  nodes  and  200  branches.  Compilation  of  the  calculating 
scheme  for  problems  of  optimizing  the  development  of  networks  is  a 
labor-intensive  process  and  I’equires  careful  preparation  of  the  input 
information  even  in  the  presence  of  check  systems.  But  the  labor 
expenditures  to  create  the  basic  information  about  the  calculating  plan 
for  optimization  calculations  of  a  given  netw^ork  are  single-  stage. 

This  information  is  used  repeatedly  for  a  long  period  in  studies  of  the 
development  of  EES.  The  calculating  plan  is  only  corrected  during  each 
study  cycle. 

Compilation  of  the  task  for  optimization  (determination  of  the  area  of 
search  for  the  optimal  solution)  and  display  of  it  in  the  calculating 
plan  of  the  network  x^redetermine  to  a  significant  degree  the  process  of 
optimization.  The  erudition  of  the  investigator  is  of  important 
significance  for  compilation  of  the  optimization  task.  A  set  of 
alternative  measures,  which  are  described  by  means  of  variables,  is 
determined  in  compilation  of  the  optimization  task.  The  problems  to  be 


solved  are  those  with  discrete  variables.  Only  two  values  are  possible 
for  each  variable:  0 — a  measure  corresponding  to  the  variable  is  not 
implemented,  and  1 — the  measure  is  implemented  according  to  the 
variable.  Each  \'ariable  is  characterized  by  an  ordinal  number  and  by  a 
set  of  branches  of  the  network,  for  vliich  the  number  of  elements  of  each 
branch  per  unit  is  increased  on  this  and  on  subsequent  steps  of  the 
optimization  process  after  implementation  of  the  variable  at  each 
step  t . 

When  postulating  the  optimization  problem,  to  be  solved  by  using  the 
QRS  E,  an  alternative  measure  should  be  determined  such  that 
implementation  of  it  would  yield  an  energy  effect  (reduction  of  losses 
in  the  network ,  reduction  of  loads  and  so  on ) . 

The  optimization  program  system  for  development  of  an  QRS  E  uses 
variables  of  two  types.  Implementation  of  variables  of  the  first  type 
displays  the  measures  on  construction  of  electric  power  transmission 
lines  over  new  routes,  higher  voltage  lines  and  so  on.  No  more  than  one 
new  element  appears  in  the  corresponding  future  branches  (with  the 
number  of  elements  in  the  primary  scheme  equal  to  zero).  Implementation 
of  variables  of  the  second  type  shows  the  construction  of  new  parallel 
LEP  or  installation  of  additional  generators  at  electric  power  plants 
with  the  same  parameters  as  existing  plants,  which  is  related  to  an 
increase  of  the  number  of  elements  both  in  existing  and  in  future 
branches.  When  compiling  the  optimization  task  using  the  ORS  E  program 
system,  one  can  impose  restrictions  on  the  possible  period  of 
implementation  of  variables,  and  the  beginning  and  end  of  the  period  of 
implementation  of  alternative  measures  must  be  indicated. 

It  becomes  necessary  in  optimization  problems  to  take  into  account  the 
dependence  between  implementation  of  measures.  For  example,  to  replace 
the  wires  on  some  electric  power  transmission  line,  the  line  must  first 
be  constructed.  Logical  restrictions  must  be  established  for  Boolean 
variables  in  these  cases.  There  are  two  main  types  of  these 
restrictions ; 

1)  Disjunction — states  in  which  more  than  one  measure  (variable)  from  a 
given  set  is  implemented  are  not  permissible.  Disjunction  is  used  to 
simulate  mutually  exclusive  measures  on  developnent  of  an  electric 
network.  This  method  of  simulation  can  be  used,  for  example,  vdien 
selecting  the  optimal  variant  of  the  design  version  of  an  electric  power 
transmission  line  (type  of  current,  rated  voltage,  mark  of  wires). 

2)  Limitations  on  the  sequence  of  implementing  measures  to  develop  the 
network.  If  j  after  i  is  given,  implementation  of  the  j-th  variable  is 
not  permissible  if  the  i-th  variable  has  not  been  implemented. 
Limitations  of  type  *’j  after  i"  are  used  to  simulate  measures  which  can 
be  implemented  only  in  a  specific  sequence.  For  example,  if  temporary 
use  of  an  electric  power  transmission  line  at  lower  voltage  with 
subsequent  conversion  of  it  to  rated  voltage  is  considered,  this 
conversion  can  be  implemented  only  after  implementation  of  the  first 
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measure — introduction  of  the  electric  power  transmission  line  at  the 
lower  voltage  compared  to  the  rated  voltage. 


The  use  of  logical  restrictions  reduces  the  dimensionality  of  the 
dynamic  problems.  Thus,  for  example,  the  presence  of  10  logical 
restrictions  with  20  variables  causes  a  reduction  of  the  number  of 
possible  states  by  approximately  a  factor  of  seven.  Compilation  of  the 
optimization  task  is  made  difficult  by  the  circumstance  that  the  number 
of  variables  in  dynamic  optimization  models  is  limited:  the  maximum 
permissible  number  is  64  variables.  Accordingly,  some  aggregation 
(consolidation)  of  measures  is  possible. 

Let  us  consider  below  the  methods  of  simulation  and  optimization  which 
are  used  in  the  ORS  E  program  system. 
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— number  of  nodes  in  calculating  plan; 

— total  number  of  branches  in  calculating  plan; 

— total  number  of  branches  that  show  electric  power 

transmission  line  (LEP); 

— total  number  of  branches  that  show  electric  power 

plants  { ES ) ; 

— ordinal  number; 

— number  of  parallel  elements  of  i-th  branch  on  t-th 
development  step, 

for  LEP — number  of  parallel  electric  power 
transmission  lines, 
for  ES — number  of  units  of  same  type; 

— number  of  calculated  sets  of  loads,  considered  on 
t-th  development  step 

— length  of  electric  power  transmission  line  for  i-th 
branch,  km; 

— flow  of  active  power  in  i-th  branch  on  step  t  for 
i-th  calculated  load,  MW; 

— coefficient  proprotional  to  effective  resistance, 
for  determination  of  acti^^e  power  losses  in  one 
element  of  i-th  branch,  1/MW; 

— coefficient  proportional  to  i-eactance,  for 
determination  of  reactive  power  losses,  1/MVA; 

— coefficient  proportional  to  conductivity,  for 

detemination  of  flow  of  active  power  in  one 
element  of  i-th  branch  (economical  conductivity  of 
branch  1 ) ,  MW; 

— longitudinal  reactance  of  generator,  percent; 

— specific  active  power  losses  to  corona  in  i-th 

branch,  kW/km; 

— power  coefficient  of  power  flux,  relative  units; 
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— specific  susceptance  of  electric  power  transmission 
line,  10®  cra/km; 

— rated  voltage  in  i-th  branch,  kV; 

— capital  investments  in  LEP,  shown  by  i-th  branch, 

dependent  on  its  length,  thousand  rubles/km; 

— capital  investments  in  i-th  branch,  independent  of 

its  length  (including  those  for  substations), 
thousand  rubles; 

— capital  investments  in  i-th  branch,  shown  by 

electric  power  plant,  independent  of  ntnnber  of 
units ,  thousand  rubles ; 

— capital  investments  for  installation  of  one  unit  in 

i-th  branch,  representing  power  plant,  thousand 
rubles ; 

— capacity  of  unit,  MW; 

— specific  fuel  expenditures,  kopecks/kW-hr; 

— number  of  hours  of  utilization  of  installed  power  of 
i-th  power  plant; 

— specific  capital  investments  for  installation  of 

reactors,  thousand  rubles /MVA; 

— specific  capital  investments  in  pressurizer  on  one 

electric  power  transmission  line,  thousand 
rubles/MVA; 

— degree  of  compensation  of  charging  power  of  power 

transmission  lines,  relative  units; 

— efficiency  of  capacitive  power  of  electric 

transmission  line,  relative  units; 

— relative  costs,  dependent  on  capital  investments 

K^^(i)  for  electric  power  transmission  lines, 

relative  units; 

— relative  costs,  dependent  on  capital  investments 

Kj^^(i)  for  power  plant,  relative  units; 

— relative  costs  for  depreciation  and  maintenance  for 

electric  power  transmission  lines,  dependent  on 
K^(i),  relative  units; 

— relative  costs,  dependent  on  capital  investments 
g^(i),  relative  units; 


2T 


Kpeak' 


3'j^(i,  t) 


3\{i,  t) 


Tr(jl 

^h6 


kk 

PH(i.  j) 


Rt 

L(j) 


A  and  k 


H 


E 


Hn 


3^a» 


n 


j?(t) 

s(t) 


e(t) 


— coefficient  for  determination  of  expenditures  for 
ccHnpensation  of  load  reactive  power  losses,  thousand 
rubles/year • MVA ; 

— relative  expenditures  for  losses  in  i-th  branch, 

dependent  on  loads,  kopecks/kW*hr; 

— relative  expenditures  for  no-load  losses  in  i-th 

branch,  kopecks/kW-hr; 

— length  of  j-th  mode,  hr; 

— dynamic  reduced  expenditures; 

— additional  criterion  for  estimation  of  disruption 

of  power  balance  conditions; 

— axlditional  criterion  for  estimation  of  violation  of 

restrictions  on  power  flows; 

— number  of  parts  of  network  unrelated  to  each  other; 
— power  unbalance  in  .i-th  part  of  network  tor  i-th 

set  of  calculated  loads,  ^]W; 

— hypothetical  resistance; 

— set  of  branches  representing  LEP  for  which 

restriction  on  power  fluxes  is  disrupted,  at  J— th 
set  of  calculated  loads; 

— capacity  per  circuit,  MW; 

— coefficients  for  calculation  of  additional  criterion 
for  estimation  of  violation  of  restriction  on  power 
fluxes  are  given  in  input  information; 

— number  of  years  during  which  capital  investments  are 
made  or  annual  costs  var.v; 

— normative  factor  of  comparative  effectiveness  of 

capital  investments; 

— normative  factor  of  reduction  of  expenditures  at 
different  times; 

— reduced  expenditures  throughout  system  during  year 

— reduced  expenditures  throughout  system  at  normal 
operating  step; 

— year  of  reduction  of  expenditures  at  different  times 

(the  first  year  of  the  calculating  period  is  used  as 
the  year  of  reduction  in  the  ORS  E  program  system) ; 
— first  year  of  t-th  fJiase; 

— coefficient  of  specific  weight  of  exTJenditures 
3  (^(t))  on  t-th  development  step  in  specific 

c 

function; 

— state  of  network  on  step  t  (combination  of 

implemented  measures  on  developnent  of  electric 
power  system  during  time  up  to  t,  inclusively); 
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f(t,  e(t))  — minimum  value  of  functional — specific  function, 

corresponding!  to  optimal  process  of  develojxnent 
during  t  steps  from  initial  state  to  some  state 


e(t) ; 


g(t,  e{t)) 
e(t,  m) 


— component  of  functional,  dependent  only  on  step  t 
and  state  e { t ) ,  but  independent  of  path  of 
transition  to  given  state  e{t); 

— states  of  development  on  step  t  with  m  implemented 
measures ; 


S  — constant,  called  critical  gradient  value; 

grad*f(t,  e(t,  m) ) — relative  maximum  increase  of  functional  for  state 

e{t,  m)  according  to  comparison  to  functional  for 
state  e(t,  m  -  1)  (relative  gradient). 


Indicators  and  Criteria  of  Variants 

Active  power  losses  AP(t,  j)  on  step  t  with  the  j-th  set  of  calculated 
loads  are  determined  as  the  sum  of  load  losses  and  no-load  losses: 


hP{t,  /)  =^Pn{t,  j)  +APxx(^,  /■). 


(1) 


The  active  power  load  losses  in  MW  are  : 


n 


APh 


PlHI,  t,j)R{i) 
c{i,  t) 


(2) 


The  active  no-load  power  losses  in  MW  are: 


"jian 

APxx(/,/)=  PK(j)-10-3.c(i,/)./(t). 


(3) 


The  reactive  power  losses  in  the  electric  power  transmission  lines  in 
MVA  are  : 


jisn 


AQ  {t,  /■)  =  E  (  ~ (0  •  l0-«  •  C  (t.  /)  •  / (i) )  . 


(4) 


The  annual  costs  are  determined  as  the  sum  of  the  costs  for  depreciation 
and  maintenance,  for  energy  losses  and  for  compensation  of  the  reactive 
power  in  thousand  rubles/year: 
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H(/)=Hk(0+h„(/)+Hq(/). 


(5) 


The  annual  costs  for  depreciation  and  maintenance  in  thousand 
rubles /year  are: 


"jian  "ac 

Hk(/)=  Hk,  ji3n (if  0  +  3c(t,  0- 

1=1  2=1 


The  annual  costs  for  depreciation  and  maintenance  for  branch  i, 
representing  the  LEP,  in  thousand  rubles/year  are: 


Hk,  Jian  (f)  t)  —  Kka  (i)  (c  (i,  /)  c  (t,  0) )  (XKa~h  {Kn  (i)  •  cc^^  + 
+'g'p(t)  ■Ua^(i)  ■  10-®(/Cp-Yp  — /CcYc)aKa)/(i)  (c{i,  t)  — c(j,  0)). 


The  annual  costs  for  depreciation  and  maintenance  for  branch  i, 
representing  the  LEP,  in  thousand  rubles/year,  are: 


HK,3c(i.  t)  =/C6ji(0  •  (c(t,  t)  -c(i,  0))  •aoji+ 

+  (^3C  (0  +  ^Ka  (0  )  •  anoCT. 


(8) 


The  annual  costs  for  losses  in  thousand  rubles/year  are: 


=  1  j=l 


(9) 


‘;i3n 


+  X  P40-3"r(LO-/(0-8760-10-Mc(U)-c(LO)). 


i=l 


The  annual  costs  for  compensation  of  reactive  power  in  thousand 
rubles/ year  are : 


”ji3n 


Hg 


i=l 


Kc  •  OiKa-'^  (i) 
COS*  (pc  (i,  t) 


maxPx,*(t,  t,  /). 
j=I,r 
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Capital  investments  are  determined  as  the  sum  of  capital  investments  in 
the  electric  power  transmission  line  and  the  power  plant: 


"jisn  "ac 

^  ^ji3n(i, /)+  t)- 


(10) 


Capital  investments  for  branch  i,  representing  the  LEP,  in  thousand 
rubles,  are: 


()  —  {Kua (i)  +  +gp{i)  •  I0~®X 

X  (/CpYp  -  A^cYc)  )•/(£)}•  (c  (/,  0  -  c  (t  0) ) .  <  1 1 ) 


Capital  investments  for  branch  i,  representing  the  power  plant,  in 
thousand  rubles ,  are : 


^3C  (t  0  ~  ■^3C  "t  ^6jr  (0  ■  (^  (t.  0  ~  ^  (t  0)  )  • 


(12) 


Dynamic  expenditures  are  determined  by  the  expression  in  thousand 
rubles/year : 


3„= 


E 


{Ea-Kz{t)+Vl{t))-S{t). 


(13) 


The  dynamic  reduced  expenditures  for  the  step  process  (the  calculating 
period  is  divided  into  t  normal  operating  steps,  i.e.,  the  loads  for  all 
calendar  years  of  the  step  t  are  equal,  while  the  scheme  and  parameters 
of  the  system  do  not  vary  during  the  step)  are  determined  by  the 
expression : 


T 

3„=  ^  3c  (« (/) )  S  (0  +3,;  (0  I )  ( I  +£„n) 
(=1 


(14) 


where  the  coefficient  of  the  specific  weight  of  expenditures  on  the  t-th 
development  step  is  determined  by  the  formula: 


S(/)=  (l+£™)«n-W-  (l+£,„)«n-«('+‘).  (  15) 


31 


The  ORS  E  program  system  permits  the  investigator  to  assign  the  value  of 
coefficients  S(t)  for  t  =  1,  2,  ....  T  according  to  his  discretion.  If 
the  values  of  S(t)  are  not  assigned,  they  are  determined  by  formula 
(15).  An  additional  criterion  of  estimating  violations  of  the  condition 
of  power  balance  on  step  t  are  determined  by  the  expression  in  thousand 
rubles/year : 


ip(i)  hk 

Z.o(f)-=  EL 

2  =  1  j=l 


Pa^ihi)  •Rl'Ac^-Tri}). 


(16) 


The  additional  criterion  of  estimating  violation  of  restrictions  for 
power  flows  on  step  t  is  determined  by  the  expression  in  thousand 
rubles/year: 


ril)  Hi) 

Zu(t)  =  EE  (.Eh  •  Kjisn  0  ^ 

j— 1 

A  /  absPLji,  t,  j)  _  j  \'' 

^  r(t)  '  c{i,  t)  ■  P^ou(i) 


(17) 


Model  of  Power  Floi-7s 

The  flow  distribution  of  active  power  in  the  ORS  E  program  system  is 
determined  by  the  minimum  power  losses  in  the  network  with  regard  to 
Kirchhoff’s  first  law  as  supplementary  conditions  (1).  The  exTiression 
for  optimal  power  flows  has  the  follow'ing  form. 


PL{l)  =  Y3{t){83{i)-b3(i))<  (18) 


where  S  (i)  and  S  (.j)  are  the  economic  potentials  for  the  initial  and 
3  3 

final  node  of  branch  1,  respectively.  The  economic  potentials  are 
determined  by  a  system  of  q  equations: 


Yg  Aa=P, 


(19) 


where  A  is  the  vector  column  of  the  nodal  economic  potentials,  P  is  the 
9 

vector  column  of  the  loads  of  the  nodes,  and  is  the  quadratic  matrix 


of  economic  conductivity. 


The  active  power  losses  in  transformers  in  system-forming  networks  are 
relatively  low  compared  to  losses  in  electric  power  transmission  lines. 
Therefore,  the  calculated  layouts  of  the  network  do  not  contain 
transformer  branches  vdien  determining  the  optimum  flow  distribution  (see 
figure).  Calculation  of  the  optimal  flow  distribution  permits 
convenient  simulation  of  parallel  operation  of  an  AC  network  with  DC 
elements  f 1 ] . 
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Calculated  Layout  of  Optimization  of  Development  of 
Energy  Association; 

1 — existing  electric  power  transmission  lines  (LEP)  of  500  kV; 

2 — existing  LEP  of  500  kV,  temporarily  using  220  kV;  3 — existing 
branches  with  several  elements;  4 — future  and  alternative  LEP  of 
500  kV;  5 — future  and  alternative  LEP  of  1,150  kV;  6 — future  and 
alternative  DC  LEP  of  1,600  kV;  7 — future  and  alternative  power 

plants 


KEY: 

1.  Regional  power  plant  2.  Existing  elements 

Target  Function 

The  following  target  function  is  used  when  searching  for  the  optimal 
variant  of  development; 


33 


(20) 


The  supplementary  criterion  of  estimating  the  violations  during  the 
entire  calculating  period  is  determined  by  the  expression: 


(21) 


The  (®S  E  program  system  takes  into  account  restrictions  on  capital 
investments  K^(t)  according  to  development  steps  t=l,  2,  ...,T.  If 

the  values  of  K_(t)  are  not  given,  it  is  assumed  that  K  (t)  =  lO^  rubles 
A  ^ 

for  t  =  1,  2,  T. 

Method  of  Optimization 

All  steps  t  =  1,  2,  ....  T  are  considered  sequentially  during 
optimization*  The  optimal  development  during  t  steps  to  state  e(t)  is 
determined  at  each  optimization  step  [2].  To  do  this,  the  follouing 
recursion  relation  is  used; 


f{t,e(t))=g(i,e(t))+  min 


«(-i)<3  (<•((-!)) 


(22) 


Only  optimal  initial  states,  which  are  local  optimums  according  to  the 
following  condition,  are  sought  during  the  next  optimization  step  t 


f(i,e{i))=  min  (23) 

.'•(Oef 


The  optimal  initial  states  (OIS)  and  the  corresponding  values  f(t,  e(t)) 
are  stored,  idiich  also  permits  one  to  determine  the  optimal  development 
on  step  t  +  1  during  t  +  1  steps  to  any  state  e(t  +  1).  The  optimal 
initial  states  on  the  next  step  t  is  sought  by  sequential  consideration 
of  the  subset  of  states  {e( t,  m) } ,  m  =  1 ,  2,  ...»  where  m  is  the  number 
of  implemented  measures. 

Because  of  the  large  dimensions  of  the  problems  of  optimizing  the  ORS  E 
program  system,  the  approximate  method  of  finding  the  optimal  initial 
state  is  used. 


3k 


In  this  case,  the  optimal  initial  states  are  selected  by  the  feature 


~grad*fit,e{t,m))^6. 


(24) 


The  relative  gradient  is  determined  by  the  formula: 


grad*/(/,  e(t,  ni))  = 


f{t,e{t,m))—  min  f(t,e(/,m-l)) 

c(i,  rn~l)^e(t,  m) 


•100%. 


(25) 


The  search  for  the  optimal  initial  states  is  stopped  if  no  optimal 
initial  states  are  found  at  given  m.  One  could  find  all  the  local 
mini mums  of  the  target  relation  at  ^  =  0  and  upon  retention  of  all 
optimal  initial  states  at  each  optimization  step  and  one  could  find  the 
global  optimum.  However,  this  is  unrealistic  due  to  the  large 
dimensions  of  the  dynamic  problems.  The  real  dynamic  optimization 
problems  can  be  solved  if  the  number  of  optimal  initial  states  (S  >  0) 
is  retained  at  each  development  step  by  retaining  only  the  optimal 
initial  states  with  the  best  value  of  the  functional  for  different 
competitive  strategies.  To  reduce  the  set  of  optimal  initial  states  in 
the  program,  one  can  also  use  the  restriction  form  above  FAA(t)  for  each 
step  t=l,  2,  ...,T.  Only  those  optimal  initial  states  e'(t)  for 
vdiich  the  following  condition  is  fulfilled  are  retained  during 
optimization 


minf(t,e(t)) 

{<•(()) 


(26) 


If  FAA(t)  is  not  given,  it  is  assumed  that  FAA(t)  =  1.4,  since 
experimental  and  theoretical  studies  show  that,  if  the  values  of  the 
functional  of  some  state  differ  from  the  minimum  value  more  than  1.4,  an 
optimal  strategy  may  not  pass  through  this  state. 

Pseudo-dynamic  optimization  is  realized  at  FAA(t)  =  1.  Thus,  the  values 
of  FAA(t)  can  be  given  in  the  range  1-1.4.  The  calculating  time 
decreases  as  the  value  of  FAA(t)  decreases,  but  the  probability  of  a 
loss  of  optimal  strategy  increases.  The  values  of  FAA(t)  for  subsequent 
optimization  steps  can  be  assigned  less  than  for  the  initial  values. 

The  maximum  possible  number  of  retained  optimal  states  is  WM  =  40.  A 
decrease  of  WM  reduces  the  time  for  optimization,  but  the  probability  of 
a  loss  of  the  optimal  variant  according  to  the  integral  target  relation 
during  the  calculating  period  increases.  Pseudo-dynamic  optimization  is 
realized  at  WM  =  1.  The  procedure  of  optimization  (22)  for  steps  t  =  1, 
2,  . . . ,  T  is  performed  repeatedly  (by  the  iteration  method)  at  #  >  0. 
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Table  1.  Optimized  Measures  (Variables) 


1 

(1)  .  ^  2 

HaHMeHOoaiiHe  MeponpnflTHft 

)  Mho/kcctbo  BeTBefi  | 

(iprpaHHMeHiiH 

Wi' 

Wr 

{ro},. 

in 

1 

2 

3* 

4* 

5* 

6* 

7* 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 
11 
12 

13 

14 

15 

16 

17 

18 

19 

20 
21 


Beoa  3neKTpocTaHmiH  b  ys.ie  14 
Beoii  3.ieKTpocTaHmin  b  yane  15 
Bbo;i  3.aeKTpocTaHLiHH  b  ya.ie  3 
Bbo;i  /ISn  kB  12 — 14 
Bboa  JlSn  500  kB  12—13 


(ii 


Bboa  ;iBy^uennoH  /T 30  500  kB  12 


(  13) 
(  14) 


31 

32 

33 
20 
21 
22 

23  . 

24,  26 
24,  25 
19 
18 


Bsoii  J13n  500  kB  3 — 12 
Bboa  500  kB  10—11  —  12 
Bboa  Jisn  500  kB  6—11—12 
Bboa  J\3Tl  500  kB  13—14^ 

Bboa  J13n  500  kB  13 — 15  \ 

Bboa  3-h  napaAAeAbHoft  J1 311  500  kB 
Bboa  4>ri  napaAAeAbHoii  J13n  500  kB  1—15  1^ 
Bboa  5-h  napaAAeAbHofi  .1311  500  kB  1— 15  ^  ^ 
Bboa  3-h  napaAAeAbHoft  J13n  500  kBUo} 

j _ 2 _ 3 

Bboa  4-ft  napaAAeAbHOil  J13n  500  kB 

j _ .2 _ 3 

Bboa  5-h  napaAAeAbHOH  .1311  500  kB  (20) 
j _ 2 _ 3 

Bboa  2-h  napa.iACAbnoH  «13n  500  kB(21) 

0  ^ _ g 

Bboa  /13n  1150  kB  6—16  (22) 

Bboa  /13n  nocTOfliiiioro  lOKa  6—16  (  2  3) 
nepeBOA  J13n  c  220  iia  500  kB  (24) 


^2—2 

2—2 

2 _ 2 

1— .3 
J-3 
1—3 
1—3 

1— 3 
a-3 

2— 3 
2—3 
1-h3 
1-3 
1—3 


12 

13 


5.  6 

— 

1-^3 

— 

5,  6 

— 

1-^3 

15 

__ 

5,  6 

— 

1—3 

16 

V- 

9,  10 

— 

1—3 

— 

28 

— 

— 

1—3 

— 

29 

_ 

— 

a-3 

— 

27 

— 

17 

1—3 

— 

Notes:  i — number  of  variable;  {h]i^ — set  of  future 
branches  for  which  the  number  of  elements  assumes  the  ^ 
value  of  1  after  implementation  of  the  i-th  measure;  {Lli^ 
set  of  existing  and  future  branches  for  which  the  number 
of  elements  increases  by  1  after  implementation  of  the 
i-th  measures;  ILoli — set  of  branches  for  which  conductivity 
is  set  equal  to  zero  (the  branch  is  disconnected)  after 
implementation  of  the  i— th  measure;  period  oi 

permissible  implementation  of  i-th  measiire  (t  is  the  ordinal 
number  of  the  optimization  step);  in —number  of  measure, 
after  implementation  of  w^hich  implementation  of  measure  i 
is  permitted. 


KEY 


1.  Name  of  measures 

2.  Set  of  branches 

3.  Restrictions 

4.  Introduction  of  power  plant  in  node  14 

5.  Introduction  of  power  plant  in  node  15 

6.  Introduction  of  power  plant  in  node  3 


(Key  continued  on  follow’ing  page) 
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7.  Introduction  of  LEP  12-14  of  500  kV 

8.  Introduction  of  LEP  12-13  of  500  kV 

9.  Introduction  of  two-circuit  LEP  of  12-13  of  500  kV 

10.  Introduction  of  LEP  3-12  of  500  kV 

11.  Introduction  of  LEP  10-11-12  of  500  kV 

12.  Introduction  of  LEP  6-11-12  of  500  kV 

13.  Introduction  of  LEP  13-14  of  500  kV 

14.  Introduction  of  LEP  13-15  of  500  kV 

15.  Introduction  of  third  parallel  LEP  1-16  of  500  kV 

16.  Introduction  of  fourth  parallel  LEP  1-16  of  500  kV 

17.  Introduction  of  fifth  parallel  LEP  1-16  of  500  kV 

18.  Introduction  of  third  parallel  LEP  1-2-3  of  500  kV 

19.  Introduction  of  fourth  parallel  LEP  1-2-3  of  500  kV 

20.  Introduction  of  fifth  parallel  LEP  1-2-3  of  500  kV 

21.  Introduction  of  second  parallel  LEP  3-5-6  of  500  kV 

22.  Introduction  of  LEP  6-16  of  1,150  kV 

23.  Introduction  of  DC  LEP  6-16 

24.  Switch  of  LEP  from  220  to  500  kV 


The  optimization  is  performed  on  the  first  iteration  at  some  given  value 
considers  all  variables.  The  value  of  8  decreases  on 

subsequent  iterations  at  given  step  A#  (S.8  <  Preoptimization  of 

na  H 

the  optimal  variant  of  development  is  performed.  The  ORS  E  program 
system  presents  the  investigator  the  opportunity  to  perform 
preoptimization  of  some  other  variant  as  well  if  desired.  In  this  case, 
the  investigator  should  indicate  the  ordinal  number  of  the  preoptimized 
variant.  If  the  ordinal  number  is  not  given,  the  optimal  variant  is 
preoptimized . 

Thus,  the  optimization  problem  is  solved  with  a  reduced  number  of 
variables  and  with  increased  accuracy  of  finding  optimal  initial  states. 
Optimization  ends  when  the  values  of  b  reach  zero,  i.e.,  when  finding 
the  exact  values  of  the  local  minimums  of  the  target  relation, 
corresponding  to  the  sought  optimal  initial  states,  it  is  guaranteed 
during  searching  for  the  optimal  initial  states. 

Example 

The  calculating  scheme  for  optimization  of  development  of  a  system- 
forming  network  with  refinement  of  the  location  of  new  power  plants  is 
presented  in  the  figure.  The  considered  energy  association  is  separated 
into  two  regions — ^RES  1  and  RES  2,  having  different  cost  of  energy.  The 
period  to  be  studied  is  divided  into  three  development  steps  with  length 
of  5  years.  A  total  of  10  calculated  sets  of  loads  is  considered  on 
each  development  step.  Introduction  of  pxDwer  plant  21  is  provided 
during  the  first  optimization  step,  and  the  problem  of  generation  of 
power  by  this  station  must  be  solved.  Another  power  plant  must  be 
introduced  on  the  second  step,  and  its  location  must  be  refined  (node 
14,  15  or  3).  The  main  problem  to  be  studied  is  that  of  selecting  the 
variant  of  electric  power  transmission  from  node  6  to  node  16.  There 
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are  three  alternatives:  1)  enlargement  of  the  500-kV  network,  2) 
construction  of  a  1,150-kV  LEP  and  3)  construction  of  a  DC  LEP.  A 
special  problem  which  is  also  studied  in  this  problem  is  determination 
of  the  moment  of  transition  of  LEP  8-17 ,  temporarily  utilized  on  220  and 
500  kV  (construction  of  a  500/220-kV  substation  at  node  8). 

Table  2.  Logic  Restrictions  of  Type  "Or" 

Set  of  variables  from  which 
implementation  of  no  more  than 
one  variable  is  permissible  in 
Number  of  restriction  one  development  variant 

1  1,  2,  3 

2  5,  6 

3  19,  20 

Table  3,  Results  of  Optimization 
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A  list  of  measures  to  be  optimized  according  to  development  of  the 
network  is  presented  in  Table  1.  Variables  12-18  are  those  of  the 
second  type.  They  are  described  by  a  set  of  existing  branches,  for 
which  the  number  of  elements  is  increased  by  one  by  implementation  of 
the  measxire.  Variables  15-18  are  described  by  two  branches,  while  the 
remaining  variables  of  second  type  are  described  by  only  one  branch. 

The  specific  sequence  of  realizing  them  is  observed  when  using  variables 
of  the  second  type  and  logic  restrictions  "after"  are  used  (Table  1, 
column  7).  Restrictions  for  the  period  of  implementation  of  alternative 
measures  (Table  1,  column  6)  and  logic  restrictions  of  the  "or'  type 
(Table  2)  are  also  used  in  the  considered  example. 

Introduction  of  these  restrictions  is  determined  by  postulation  of  the 
optimization  problem.  For  example,  those  layouts  of  the  network  are 
impossible  in  which  branches  21  (variable  5)  and  22  (x'ariable  6)  are 
present  simultaneously,  since  these  branches  represent  the  same  line  in 
a  one-  or  two-circuit  version.  Implementation  of  more  than  one  measure 
from  set  1,  2,  3  is  also  impossible,  since  only  one  new  power  plant  is 
planned  for  construction  during  the  second  development  step. 

The  results  of  optimization  are  presented  in  Table  3. 
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